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ABSTRACT 


Data  from  a  research  tower  in  Lake  Ontario  are  used  to  study  the  validity  of  Monin- 
Obukhov  scaling  in  the  marine  atmospheric  boundary  layer  tmder  various  wave  con¬ 
ditions.  It  is  found  that  over  pure  wind  seas,  the  velocity  spectra  and  cospectra  follow 
established  universal  scaling  laws.  However  in  the  presence  of  swells  outrunning  weak 
winds,  velocity  spectra  and  co-spectra  no  longer  satisfy  universal  spectral  shapes.  Here, 
Monin-Obukhov  similarity  theory,  and  the  classical  logarithmic  boundary  layers  are 
no  longer  valid.  It  is  further  shown  that  in  the  presence  of  such  swells  the  momentum 
flux  can  be  significantly  modified  over  pure  wind  sea  values.  The  implications  of  these 
findings  for  bulk  flux  estimations  and  on  the  inertial  dissipation  method  for  calculating 
fliixes  are  discussed. 

Keywords:  Air-sea  fluxes,  wind  stress,  wind  waves,  swell,  universal  spectra,  Monin- 
Obukhov  theory 
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1,  Introduction 

The  basis  of  most  recent  work  on  turbulence  in  the  atmospheric  surface  layer  lies 
in  the  ideas  first  proposed  by  Obukhov  (1946)  and  Monin  and  Obukhov  (1954).  In 
stationary  and  homogeneous  conditions,  the  stress  r  is  assumed  to  be  constant  in  a 
region  near  the  interface,  but  well  above  the  viscous  sublayer.  Here,  the  stress  is  given 

by 

T  =  —pu'w'  =  (1) 


where  u'  and  w'  are  the  horizontal  (in  line  with  the  mean  wind)  and  vertical  compo¬ 
nents  of  the  turbulent  fiuctuations  of  the  wind  field,  p  is  the  air  density,  and  the  overbar 
represents  time  averaging.  In  Monin-Obukhov  (MO)  similarity  theory  the  friction  ve¬ 
locity,  is  introduced  as  a  velocity  scale.  In  neutral  conditions,  with  the  assumption 
of  eddy  viscosity  increasing  linearly  with  height  z  and  friction  velocity,  this  leads  to  the 
well  known  logarithmic  velocity  profile  of  a  shear  layer. 

More  generally,  buoyancy  is  an  additional  source  (or  sink)  of  turbulence.  An  important 
consequence  of  Monin-Obukhov  (MO)  similarity  theory  is  that  buoyancy  can  be  brought 
into  the  above  picture  through  the  addition  of  a  single  length  scale  L,  the  Obukhov 
length: 


L  =  * 

Kg(w'd'/To  -1-  O.Slw'q'Y 

where  k  is  the  dimensionless  von  Karm^  constant,  g  is  the  gravitational  constant, 
To  the  reference  absolute  temperature,  and  &  and  q'  are  the  turbulent  fluctuations 
of  potential  temperature  and  specific  humidity  respectively.  The  Obukhov  length 
represents  the  height  at  which  mechanical  (shear)  and  buoyant  forcing  are  equal.  It 
follows  from  MO  theory  that 


dU  KZ 

dz  u* 


=  MO 


where  C  =  zjL.  Integrating  (3),  the  velocity  profile  is 

^  =  -K— )  -  V’«(C)]- 
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where  Uz  is  the  mean  wind  speed  at  z  meters,  Us  the  mean  surface  velocity,  Zo  is  the 
roughness  length,  where  Uz^  =  Us,  and 

<!>.  =  /  IJ— (5) 

•'Co  s 

In  neutral  conditions,  C  ^  0,  <^u  — >  1  and  V’u  — ^  0,  yielding  the  classical  logarithmic 
profile  shape. 

Since  the  introduction  of  MO  theory,  considerable  work  has  been  carried  out 
to  determine  the  dimensionless  function  and  the  von  Karman  constant.  Early 
experimental  verification  of  MO  similarity  in  the  atmosphere  over  land  was  given  by 
Zilitinkevich  and  Chalikov  (1968),  Dyer  and  Hicks  (1970)  and  Businger  et  al.  (1971). 
Hogstrom  (1996)  provides  a  summary  of  more  recent  work.  If,  as  suggested  by  Yaglom 
(1977),  the  various  results  are  normalised  so  that  k  =  0.4,  Hogstrom  (1996)  finds  good 
agreement  amongst  many  experiments. 

Although  very  little  profile  data  over  water  have  been  presented  to  date,  it  is  generally 
assumed  that  the  profile  relations  over  the  ocean  are  identical  to  those  over  land  (c.f. 
Large  and  Pond  1981,  Yelland  and  Taylor  1996).  However,  it  has  long  been  recognised 
that  the  physics  at  the  air-sea  interface  are  significantly  different  from  those  at  an  air¬ 
land  interface  (Stewart  1961).  Over  the  sea,  energy  and  momentum  are  transferred  from 
the  winds  to  the  surface,  resulting  in  waves  and  surface  currents.  As  the  waves  evolve 
the  surface  roughness  changes,  hence  the  air  fiow  over  the  ocean  surface  is  coupled  with 
the  evolution  of  the  surface  itself.  In  p^ticular,  near  the  sea  surface  (again,  away  from 
the  viscous  sublayer)  the  stress  is  ma4e  up  not  only  of  a  turbulent  component,  but  also 
of  a  wave-induced  component,  r’",  where  this  latter  component  accounts  for  the  transfer 
of  momentum  to  the  waves  (Janssen  1989).  Away  from  the  surface  the  turbulent  stress 
dominates,  while  the  wave  induced  stress  is  increasingly  important  near  the  surface. 
Given  that  the  total  stress,  t  =  —  pu'w\  is  constant,  the  increase  in  r’®  near  the 

STirface  is  balanced  by  a  reduction  in  the  turbulent  stress.  Hence  in  this  region  the 
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velocity  profiles  (supported  by  the  turbulent  stress)  are  expected  to  deviate  from  their 
over-land  shapes  (4)  (Stewart  1961). 

Using  an  empirical  wave  model,  Janssen  (1989)  estimated  the  thickness  of  the 
wave  boundary  layer  (WBL,  where  the  wave  induced  stress  is  significant)  to  be 
u\lg  «O(0.1m)  for  young  waves;  less  for  mature  waves.  Makin  and  Mastenbroek  (1996) 
estimate  r"*  to  be  under  5%  of  the  total  stress  at  the  height  of  1  m,  under  1%  by  3  m 
(for  U  =  15  ms“^).  These  estimates,  made  assuming  pure  unimodal  wind  seas,  indicate 
that  the  WBL  lies  below  typical  measurement  heights  in  the  field.  Hence,  we  might 
expect  that  measured  profiles  over  the  sea  would  be  similar  to  those  over  land.  Indeed, 
recent  open  ocean  results  of  Edson  and  Fairall  (1998),  based  on  data  estimated  to 
have  been  taken  above  the  WBL,  show  strong  support  for  MO  similarity  in  the  marine 
atmospheric  surface  layer.  There  are,  however,  indications  that  this  is  not  always  the 
case.  Holland  tt  al.  (1981)  presented  \mpublished  data  from  the  Donelan  et  al.  (1974) 
experiment  showing  profiles  where  the  wind  speed  increased  from  a  height  of  12.3  m 
down  to  2.4  m.  In  the  case  presented,  relatively  light  winds  were  blowing  in  the  same 
direction  as  a  strong  swell.  Similar  results,  under  similar  conditions,  were  reported 
recently  by  Smedman  et  al.  (1999).  Hence,  departures  from  MO  similarity  have  been 
observed  over  waves. 

This  has  particular  significance  with  respect  to  efforts  to  determine  the  so-called 
bulk  relations  by  which  one  can  estimate  surface  fluxes  from  more  readily  measured 
values,  such  as  mean  wind  speeds  and  temperatures.  This  is  important  in  improving  the 
accuracy  of  atmospheric  and  oceanic  models,  which  are  often  limited  by  sparseness  of 
reliable  boundary  data.  The  drag  coefficient,  relating  the  momentum  flux  to  the  mean 
wind  speed,  is  Cdz  =  —  Us)^)-  According  to  MO  theory,  Cdz  =  f(z/zo,C), 

which  follows  from  (4).  It  is  standard  practice  to  eliminate  the  dependence  on  stability 
and  height  by  considering  the  neutral  drag  coefficient  at  a  reference  height  of  10m, 
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CdioN  =  'rl{p{UioN  -  Us)'^),  where  U^n  =  Uz  +  (u*/k)V>«(C)-  Hence,  CmoN  is  expected 
to  be  a  function  of  the  roughness  alone. 

Using  dimensional  arguments,  Charnock  (1955)  proposed  the  relation  Zo  =  ocul/g, 
where  a,  the  Charnock  constant,  is  of  order  0.01.  The  resultant  dependence  of  Cdion 
on  wind  speed  was  confirmed  by  Smith  (1980).  More  recently  it  has  been  shown  that 
Cdion  is  also  dependent  on  wave  age  (Donelan  1990;  Smith  et  al.  1992;  Donelan  et 
al.  1993),  implying  that  a  is  in  fact  a  function  of  wave  age.  This  was  first  proposed  by 
Kitaigorodskii  and  Volkov  (1965). 

However,  it  is  evident  that  scatter  in  experimental  determinations  of  Cdion  can 
not  be  explained  by  variability  in  wind  speed  and  wave  age  alone.  One  factor  which 
has  been  proposed  to  account  for  some  of  the  observed  scatter  is  swell.  Geemaert  et 
al.  (1986),  Donelan  et  al.  (1997)  and  Rieder  (1997)  all  report  increased  variability  in 
the  drag  coefficient  in  the  presence  of  background  swells.  Donelan  et  al.  (1997)  and 
Drennan  et  al.  (1999)  presented  evidence  showing  that  strong  swells  (those  whose 
energy  greatly  exceeds  that  of  the  wind  sea)  running  against  light-moderate  winds 
were  associated  with  CdioN  values  significantly  enhanced  over  pure  wind  sea  values.  In 
contrast,  Dobson  et  al.  (1994)  and  Yelland  and  Taylor  (1996)  reported  no  effect  of  swell 
on  the  drag  coefficient.  Smith  et  al.  (1996),  in  their  review  of  air-sea  fluxes,  conclude 
that  the  coupling  of  the  stress  and  swell  remains  an  area  for  future  research,  especially 
in  the  open  sea. 

Here  we  investigate  the  hypothesis  that  the  presence  bf  swell  can  significantly  alter 
the  drag  coefficient,  and  that  these  changes  cannot  be  modelled  within  the  framework  of 
traditional  MO  theory.  In  particular  in  these  conditions,  in  which  the  swell  introduces 
into  the  problem  an  additional  length  scale,  MO  theory  would  appear  to  fail.  The 
outline  of  the  paper  is  as  follows.  In  section  2  we  describe  an  experiment  in  which 
eddy-correlation  stress  and  other  supporting  data  were  collected  at  a  tower  in  Lake 
Ontario.  In  section  3,  the  analysis  methods  are  described.  In  the  following  sections. 


-6- 


resTilts  axe  presented  for  both  pure  wind  sea  and  mixed  seas  (swell  plus  wind  sea).  Some 
implications  of  the  results  are  discussed  in  section  6,  and  conclusions  are  presented  in 
section  7. 

2.  The  WAVES  experiment 

During  the  autumn  seasons  of  1985-1987,  the  WAVES  (Water- Air  Vertical  Exchange 
Study)  experiments  were  carried  out  from  the  National  Water  Research  Institute 
(NWRI)  research  tower  in  Lake  Ontario.  The  tower  was  designed  specifically  for  air-sea 
interaction  research,  with  minimal  disturbance  to  both  wind  and  waves  (see  Figture  1).  It 
is  located  in  12  m  of  water,  1100  m  from  the  western  shore  of  the  lake  (Figure  2),  and  is 
connected  to  an  onshore  trailer  via  underground  cable.  Wave  heights  at  the  tower  have 
been  measured  in  excess  of  6  m  in  strong  easterly  (long  fetch)  winds,  whereas  heights  of 
less  than  1  m  axe  typical  in  strong  westerly  (short  fetch)  blows. 

During  the  WAVES  experiments,  the  tower  was  equipped  to  study  the  turbulence 
budgets  on  both  sides  of  the  interface.  A  full  description  of  the  air-side  apparatus  is 
provided  in  Donelan  et  al.  (1999),  while  the  water-side  turbulence  measurements  are 
discussed  in  Terray  et  al.  (1996). 

We  describe  here  the  instruments  relevant  to  this  study.  A  Gill  anemometer-bivane 
(R.M.  Young)  was  mounted  on  a  mast  on  the  tower  (Figure  1),  and  yielded  the  three 
components  of  the  wind  vector,  and  thence  the  wind  stress  vector,  at  a  height  of 
approximately  12  m.  The  bivane’s  anemometer  was  calibrated  using  the  100  m  long 
towing  carriage  current  meter  calibration  facility  at  NWRI.  The  signals  from  the  bivane 
anemometer  were  corrected  for  frequency  response  of  both  the  styrofoam  propeller 
(Hicks  1972)  and  vane  (MacCready  and  Jex  1964).  Wind  tunnel  tests  carried  out  at 
NWRI  indicate  a  propeller  distance  constant  of  1  m.  The  damping  coefficient  and 
natxiral  wavelength  of  oscillation  for  the  second  order  vane  correction  were  found  to  be 
0.52  and  5.8  m  respectively.  The  separation  distance  between  the  propeller  and  vane  is 
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0-73  m.  Sensors  for  temperature  and  humidity  were  also  located  on  the  mast  at  about 
12  m. 

Near  the  surface,  high  frequency  velocity  fluctuations  were  measured  with  an  X- 
film  anemometer  (Thermal  Systems,  Inc.).  The  two  cross  channels  were  converted  to 
horizontal  and  vertical  velocity  components  using  an  iterative  routine  based  on  King’s 
law.  The  X-film,  along  with  Pitot  and  Elliott  pressure  probes,  were  mounted  on  the 
“wind  profiler”,  an  instrumented  support  structure  designed  to  move  4  m  along  a  track 
parallel  to,  and  2.5  m  outboard  of,  the  northeast  leg  of  the  tower  (41  cm  dia).  The 
profiler,  2.5  m  in  height,  was  vaned  so  that  it  would  point  into  the  wind,  and  positioned 
relative  to  the  surface  by  computer  control  from  shore.  Further  details  on  the  profiler, 
and  on  the  operation  and  calibration  of  the  X-film  and  pressure  probes,  are  available  in 
Donelan  et  al.  (1999). 

The  surface  elevation  was  measured  using  capacitance  waves  gauges.  The  directional 
distribution  of  the  wave  field  was  determined  using  maximum  likelihood  analysis  (Capon 
1969)  on  the  signals  from  an  array  of  six  gauges  forming  a  centered  pentagon  of  25  cm 
radius.  The  near  surface  current  was  measured  at  two  depths,  roughly  0.5  and  2m,  using 
drag  sphere  current  meters  (Donelan  and  Motycka  1978). 

3.  Analysis 

The  WAVES87  campaign  took  place  during  the  eight  week  period  between  late 
October  and  mid  December  1987.  During  this  period,  306  flux  runs  of  30-min  were 
gathered.  Meteorological  conditions  dining  these  runs  are  summarised  in  Figure  3. 
During  the  autumn,  the  area  is  dominated  by  three  types  of  events.  The  prevailing 
winds  are  SW,  and  winds  from  this  direction  clockwise  through  north  result  in  fetch 
limited  conditions  with  inverse  wave  ages  {U$  fcp)  typically  around  3-5.  Here  Cp  is 
the  phase  velocity  of  waves  at  the  peak  of  the  spectrum,  and  U$  the  component  of 
the  wind  in  peak  direction  of  the  waves  (Donelan  et  al.  1985).  Below,  we  denote  the 
inverse  wave  age  as  U/cp.  Significant  wave  heights  range  from  well  under  0.5  m  for  pure 
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wind  sea  conditions  to  over  2  m  in  the  case  of  counter  swells  propagating  down  the 
laJce.  Winds  from  the  east  yield  waves  of  order  1  m,  with  inverse  wave  ages  from  near 
zero  to  2.  In  these  cases,  it  is  often  difficult  to  discern  the  active  wind  sea  component 
of  the  spectrum,  especially  if  the  wind  is  decaying.  Finally,  there  are  a  significant 
nmnber  of  cases  of  NNE  winds.  In  these  cases,  the  wind  sea  is  typically  accompanied  by 
waves  propagating  down  the  lake,  from  the  direction  of  the  longest  fetch,  roughly  ENE 
(Donelan  et  al.  1985). 

The  306  runs  were  initially  examined  based  on  several  criteria.  Runs  with  mean  wind 
speeds  below  3.5  ms~^or  any  1-min  average  wind  speed  below  2  ms“  Vere  rejected  on 
the  grounds  that  the  bivane  is  not  expected  to  perform  well  in  these  conditions.  Also, 
nms  that  exhibited  large  changes  in  either  wind  speed  or  wind  direction  were  rejected  as 
being  nonstationary.  Data  were  screened  for  spikes.  These  criteria  together  eliminated  57 
nms  from  the  data  set.  In  carrying  out  the  eddy-correlation  stress  calculations,  the  time 
series  were  first  detrended,  but  no  filtering  was  carried  out. 

Over  water,  where  the  roughness  is  not  directionally  homogeneous,  the  stress  is  given 
by 

r  =  p(—u'w'i  —  v'w'j).  (6) 

Here,  u',  v'  and  w'  are  the  turbulent  components  of  horizontal  in-line  (with  the  mean 
wind),  horizontal  cross,  and  vertical  air  velocities  respectively.  For  brevity,  we  henceforth 
omit  the  primes  from  the  notation.  (6)  differs  from  (1)  in  accoimting  for  contributions  to 
the  stress  which  are  not  aligned  with  the  mean  wind  (Geemaert  et  al.  1993).  The  near 
surface  currents  were  calculated  from  velocity  data  from  the  drag  sphere  current  meters. 
For  most  runs,  the  upper  current  meter,  at  approximately  0.5  m  depth  was  used. 
However,  tmder  high  wave  conditions,  this  instrument  was  covered  to  avoid  damage; 
the  lower  current  meter  (2m)  was  used  in  these  circumstances.  For  each  run,  the  drag 
spheres  were  oriented  into  the  dominant  waves  to  yield  along- wave  and  vertical  current 
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components.  During  the  WAVES87  campaign,  along-wave  mean  cmrents  are  of  0(0.1 
ms“^) 

Drag  coefficients,  both  12m  and  lOm-neutral,  are  calculated  from  the  data.  In 
determining  the  Monin-Obukhov  length  (2),  an  iterative  routine  is  used,  u*  is  the 
measured  friction  velocity,  while  the  heat  and  moisture  fluxes  are  determined  from 
bulk  relations,  using  measured  mean  values  of  temperature  and  humidity,  and  assuming 
Dalton  and  Stanton  numbers  of  0.0012  (Large  and  Pond,  1981). 

The  runs  are  classifled  according  to  wave  conditions.  Peak  frequencies  of  both  the 
full  spectrum  and  wind  sea  (if  different)  were  found,  and  Cp  and  the  wave  age  were 
calculated.  Note  that,  in  general,  swells  observed  at  the  tower  tend  to  come  from  the 
direction  of  longest  fetch,  about  70  deg,  ENE  (Donelan  et  al.  1985).  Based  on  the 
directional  spectra,  seven  categories  were  used:  i)  pure  wind  sea,  short  fetch;  ii)  pure 
wind  sea,  long  fetch;  iii)  wind  sea  (NNE)  with  long  fetch  waves  propagating  from  the 
ENE;  iv)  short  fetch  with  dominant  counterswell;  v)  short  fetch  with  nondominant 
counterswell;  vi)  short  fetch  winds  (NW  or  S)  with  cross  swells;  and  vii)  wind  sea  (E) 
with  following  swells. 

For  comparison,  we  include  open  ocean  data  collected  from  a  small  ship  during  the 
SWADE  experiment  (Donelan  et  al.  1997).  These  are  eddy  correlation  data  from  a  K- 
Gill  anemometer  mounted  on  a  mast  near  the  bow,  with  full  motion  correction  carried 
out  using  the  measured  six  degrees  of  ship  motion.  Directional  wave  spectra  were 
obtained  from  the  ship,  using  an  array  of  wave  staffs,  again  taking  into  account  the 
motion  of  the  ship. 

4.  Pure  wind  sea  data 

In  Fig.  4  we  plot  12  m  non-neutral  drag  coefficients  versus  wind  speed  for  the  full 
data  set.  The  large  degree  of  scatter  is  typical  of  such  plots,  and  points  to  the  likelihood 
that  variables  other  than  wind  speed  are  important  in  parameterising  Cj).  These  data, 
and  all  data  discussed  below,  are  from  the  12m  bivane,  imless  otherwise  specified.  As  a 
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starting  point  to  further  analysis,  we  restrict  attention  to  cases  of  ptire  wind  sea.  This 
eliminates  possible  variability  due  to  the  presence  of  background  swells  (Geernaert  et 
al.  1986;  Donelan  et  al.  1997;  Rieder  1997).  The  wind  sea  data,  107  runs,  are  denoted 
by  circles  in  Fig.  4,  with  the  short-  and  long-fetch  data  denoted  by  hollow  and  filled 
circles  respectively.  The  scatter  is  not  significantly  reduced  from  the  full  data  set, 
but  a  strong  stratification  between  the  short-fetch,  0(1  km)  and  long-fetch,  0(100 
km),  data  is  evident.  This  is  the  wave  age  effect  reported  by  Donelan  (1990),  Smith  et 
al.  (1992)  and  others.  The  short  fetch  data  are  primarily  young  wind  sea,  mostly  with 
2.5  <  Ufcp  <  4.5,  while  the  long  fetch  data  represent  more  mature  seaa,  with  inverse 
wave  ages  around  1  to  2  (Table  1). 


vie. 

long  fetch 

short  fetch 

SWADE 

0.8  -  1.5 

27 

5 

15 

1.5  -  2.5 

11 

9 

8 

2.5  -  3.5 

0 

33 

0 

3.5  -  4.5 

0 

22 

0 

Table  1.  Pure  wind  sea  runs  classified  by  inverse  wave  age,  U Icp.  Data  are  WAVES- 
long  fetch,  WAVES-short  fetch,  and  SWADE,  Donelan  et  al.  1997. 

An  alternative  means  of  showing  this  is  via  the  roughness  length,  z©,  which  is  related 
to  the  drag  coefficient  as  discussed  above.  Calculating  the  dimensionless  surface 
roughness,  that  is  the  ratio  of  Zo  to  <7,  the  standard  deviation  of  the  wave  height,  we  see 
in  Figure  5  that  Zolcr  increases  with  inverse  wave  age  as  measured  by  Ut^/cp.  Younger 
waves  are  rougher,  and  the  drag  coefficient  over  them  is  commensurately  higher.  The 
pture  windsea  data  set  is  in  very  good  agreement  with  the  Donelan  (1990)  curve,  which 
was  established  at  the  same  site,  but  with  data  gathered  in  a  different  year,  viz.  1976. 
The  SWADE  open  ocean  data,  with  wave  ages  similar  to  the  long  fetch  data  (Table  1) 
are  also  seen  to  support  the  curve. 
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If  we  now  focus  on  the  mature  pure  wind  sea  data,  long  fetch  plus  SWADE,  we  find 
the  scatter  oi  Cd  versus  U  much  reduced  (r  =  0.73,  versus  r  =  0.41  for  all  WAVES 
data).  In  Fig.  6,  the  long-fetch  and  SWADE  wind  sea  10  m  neutral  drag  coefficients  are 
plotted.  The  two  data  sets  are  in  very  good  agreement,  and  a  regression  to  them  yields 
lOOOCuioAT  =  0.117C/’ioj\r  +  0.15,  with  a  correlation  coefficient  of  0.8.  For  comparison,  the 
curve  of  Yelland  and  Taylor  (1996)  is  plotted:  clearly  the  low  wind  speed  behaviour  of 
the  Yelland  and  Taylor  data  is  not  observed  here. 

A  consequence  of  Monin-Obukhov  similarity  theory  is  that,  in  some  wavenumber 
range,  we  expect  the  velocity  spectra  to  follow  some  universal  shape  when  scaled  in  the 
appropriate  manner  (Busch  1973).  Considerable  effort  has  been  spent  in  determining 
the  universal  forms,  and  there  are  several  alternative  formulations.  Kaimal  et  al.  (1972), 
using  data  profiles  from  the  1968  Kansas  experiment,  normalised  the  velocity  spectra 
by  spectral  energy  in  the  inertial  subrange.  Their  classic  results  indicate  strong 
support  for  universal  spectral  shapes,  with  the  dimensionless  spectra  being  functions  of 
dimensionless  frequency,  and  stability,  z/L.  The  spectra  of  Miyake  et  al.  (1970), 

based  on  near-neutral  data  over  the  coastal  sea,  demonstrate  similar  xmiversality 
over  the  ocean.  Other  marine  data  sets,  e.g.  Smith  and  Anderson  (1984),  support  the 
proposed  curves. 

In  Fig.  7,  we  plot  universal  Sww,  Sun  and  Suw  spectra,  calculated  from  the  pure  wind 
sea  data  (all  wave  ages),  using  the  scaling  of  Miyake  et  al.  (1970).  The  curves  of  Miyake 
et  al.  ,  derived  from  visual  fits  to  their  plots,  are  also  shown  (dashed).  For  Sww  and 
Suwt  the  Miyake  curves  fall  within  one  standard  error  (shaded  area)  of  the  WAVES87 
miiversal  curve  except  at  high  frequency,  where  the  Miyake  et  al.  curves  depart  from 
the  expected  inertial  subrange  behaviour  (due  to  limited  anemometer  response  —  Kaimal 
et  al.  1972).  For  the  horizontal  velocity  spectrum  Suu-,  the  WAVES87  imiversal  cmrve 
is  higher  than  that  of  Miyake  et  al.  by  nearly  a  factor  of  two  at  lower  frequencies, 
likely  due  to  different  mesoscale  behaviour  at  the  two  sites.  For  the  WAVES87  data 
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(Ttt/cTu)  =  2.8  ±  0.5  (showing  1  standard  deviation),  compared  to  1,9  ±  0.2  for  the  Miyahe 
et  al.  data. 


Run 

Time 

JD 

length 

min 

U 

n,c-l 

ms 

—uw 

2  -2 
ms 

—vw 

2  -2 
ms 

2  —2 
ms 

(Ty; 

ms 

<Ty,/U 

Ta 

°C 

Ty, 

°C 

m 

fp 

Hz 

Vjc, 

Symbol 

a)  Pure  wind  sea  runs 

87128 

331.55 

90 

9.60 

0.112 

-0.003 

0.100 

0.38 

0.039 

2.24 

4.52 

1.45 

0.18 

1.15 

0 

» 

» 

n 

8.14 

0.121 

0.073 

0.28 

0.034 

Hotfilm  at  2.1m 

* 

87176 

344.79 

90 

4.70 

0.019 

-0.003 

0.025 

0.22 

0.046 

5.44 

4.20 

0.10 

0.85 

2.47 

X 

b)  Swell  runs,  low  wind 

speed 

87056 

319.85 

90 

4.30 

0.017 

-0.002 

0.028 

0.14 

0.033 

7.03 

6.85 

0.67 

0.22 

0.61 

o 

87144 

333.33 

90 

4.76 

-0.006 

-0.004 

0.028 

0.19 

0.040 

5.00 

4.13 

1.45 

0.16 

0.48 

> 

87145 

333.41 

90 

4.92 

-0.001 

-0.004 

0.032 

0.17 

0.034 

4.83 

4.20 

1.44 

0.16 

0.44 

♦ 

87146 

333.64 

90 

3.90 

-0.012 

0.003 

0.006 

0.09 

0.024 

5.51 

4.30 

1.17 

0.16 

0.44 

• 

87147 

333.72 

60 

4.27 

-0.001 

-0.001 

0.019 

0.13 

0.030 

5.18 

4.33 

1.11 

0.16 

0.47 

0 

c)  Swell  runs,  moderate  wind  speed 

87117 

330.01 

60 

7.49 

0.018 

-0.009 

0.053 

0.26 

0.035 

2.88 

4.96 

2.53 

0.13 

0.74 

X 

87119 

330.26 

90 

7.57 

0.059 

0.039 

0.064 

0.31 

0.041 

0.09 

4.88 

1.46 

0.14 

0.60 

* 

87120 

330.39 

60 

8.01 

0.078 

0.055 

0.074 

0.35 

0.043 

0.13 

4.75 

1.27 

0.15 

0.63 

0 

Table  2:  Summary  of  WAVES87  data  highlighted  in  the  text.  U  is  the  mean  wind 
speed,  and  u,  v  and  w  are  the  turbulent  components  of  horizontal  in-line  (with  the  mean 
wind),  horizontal  cross,  and  vertical  air  velocities  respectively.  u^-ID  is  the  square  of 
the  friction  velocity  calculated  via  the  inertial  dissipation  method  of  Large  and  Pond 
1981,  and  <Trv  is  the  standard  deviation  of  w.  Ta,  Ty,  and  Ha  are  the  mean  air  and  water 
temperatures,  and  the  significant  wave  height,  respectively,  fp  and  Cp  are  the  frequency 
of  the  peak  of  the  surface  wave  spectrum,  and  the  corresponding  phase  velocity.  Unless 
otherwise  specified,  data  are  from  a  bivane  at  12m.  The  symbols  in  parts  (a),  (b)  and 
(c)  refer  to  Figures  8,  11-14,  and  15  respectively. 
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The  data  axe  also  categorised  by  stability,  and  the  spectra  in  four  stability  ranges  are 
shown.  There  is  little  evidence  for  a  departure  from  the  singe  universal  curve  based 
on  stability.  Only  the  most  stable  data  (circles)  show  departure  from  the  imiversal 
behaviour,  but  with  only  4  runs,  the  significance  is  low.  In  Fig.  8,  we  plot  typical 
spectra  from  two  90-min  runs:  one  with  light-moderate  12m  winds  (4.7  ms”^)  and  short- 
fetch  waves  {Ujcp  =  2.47;  =  0.1m);  the  other  with  stronger  12m  winds  (9.6  ms~^) 

and  long-fetch  waves  {UJcp  =  1.15;  Hs  =  1.4m)  -  see  Table  2a.  The  plot  shows  scatter 
arotmd  the  imiversal  curve  typical  of  the  pure  wind  sea  data:  in  general  there  is  very 
good  agreement.  In  addition  to  the  12m  bivane  data,  data  from  the  X-film  anemometer 
at  2.1  m  are  shown  for  the  long  fetch  case  {Hs  =  1.4  m).  Again,  these  data  are  seen  to 
support  the  universal  curve.  This  is  consistent  with  the  estimates  of  Janssen  (1989)  and 
Maldn  and  Mastenbroek  (1996)  as  discussed  above,  and  shows  strong  support  for  the 
validity  of  MO  similarity  over  pure  wind  seas  over  a  wide  range  of  wave  age. 

5.  Stress  data:  Mixed  wind  sea  and  swell 

The  stress  angle  has  been  shown  to  deviate  significantly  from  the  wind  direction  in 
some  situations,  and  this  has  been  attributed,  at  least  in  part,  to  the  effect  of  swells 
(Geernaert  et  al.  1993;  Rieder  et  al.  1994).  In  Figure  9,  we  plot  the  relative  angle  of  the 
stress  vector  to  the  mean  wind  (u)  direction,  arctan(tni;/t{W).  When  swell  is  present, 
particularly  at  low  wind  speeds,  there  is  often  a  considerable  difference  between  the 
wind  and  stress  directions.  For  the  pure  wind  sea  data,  the  mean  stress  direction  is 
generally  in  line  with  the  wind. 

In  studying  the  effects  of  swell  on  the  stress,  we  focus  on  the  long-fetch  data  so  as 
to  remove  the  additional  effect  of  wave  age.  In  Fig.  10  the  long-fetch  WAVES  data  axe 
plotted,  classified  as  either  pure  wind  sea  or  wind  sea  with  following  swell.  For  contrast, 
a  subset  of  the  SWADE  data  (Donelan  et  al.  1997)  are  also  plotted:  the  pure  wind 
sea  cases  and  the  cases  of  wind  sea  with  counter  swell.  Whereas  the  SWADE  counter- 
swell  runs  were  observed  to  result  in  enhanced  stress,  the  following  swell  cases  here  are 
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seen  to  be  associated  with  a  distinct  decrease  in  stress.  These  data  include  cases  where 
the  stress  is  near-zero,  and  sometimes  negative  (i.e  upward).  Hence,  the  WAVES  data 
provide  strong  support  for  a  dependence  of  stress  on  swell  properties. 

In  Figure  11  horizontal  and  vertical  velocity  spectra  from  the  bivane,  and  the  surface 
wave  height  spectrum  are  plotted  for  a  typical  diminished  or  upward  flux  case  (run 
87146  -  see  Table  2b).  These  are  runs  with  ‘fast’  swells,  Ujcp  »  0.45,  travelling  in 
roughly  the  same  direction  as  the  wind.  At  the  bivane  height  of  12  m,  the  signatmre  of 
the  swell  {Hs  =  1.2m)  is  clearly  evident  in  the  wind  velocity  spectra. 

In  Figure  12  are  plotted  spectra  of  bivane  vertical  and  horizontal  (in-line)  velocity 
for  following  swell  cases,  along  with  the  universal  curve  derived  from  the  pure  wind  sea 
data.  For  the  purposes  of  comparison  on  the  dimensional  plot,  the  data  are  restricted 
to  following  swell  in  the  mean  wind  speed  range,  U  »  3.9  —  4.9ms“^.  The  data  are 
outlined  in  Table  2b.  Typically,  U jcp  w  0.45  and  ITs-swell  >•  ifs-windsea:  these  are 
examples  of  fast,  strong  swell.  As  is  expected,  the  data  stratify  to  some  extent  by  wind 
speed:  the  filled  and  hollow  circles  represent  the  lowest  and  next-lowest  wind  speed  runs 
respectively.  The  universal  curve  is  rendered  dimensional  using  the  values  from  a  piire 
wind  sea  run  in  the  same  wind  speed  range  (run  87176,  Table  2a).  In  Fig.  12  it  is  clear 
that  while  turbulence  levels  in  the  inertial  subrange  follow  the  pure  wind  sea  data  quite 
well,  there  is  considerable  departure  at  lower  frequencies,  most  notably  at  frequencies 
lower  than  that  of  the  swell  waves  where  the  turbulence  levels  are  reduced  by  up  to  an 
order  of  magnitude  over  pure  wind  sea  values. 

In  Figure  13,  the  same  data  are  plotted,  this  time  in  the  universal  coordinates  of 
Miyake  et  al.  (1970).  Here,  ctu,  is  determined  from  the  turbulent  velocity  components 
alone,  i.e.  excluding  the  wave-coherent  component.  It  is  evident  that  spectral  similarity 
does  not  hold  when  the  swells  are  present.  The  inertial  subrange  values  are  consistently 
higher  than  those  of  the  pure  wind  sea  cases,  while  the  pure  wind  sea  ciurves  overesti¬ 
mate  the  energy  at  low  frequencies.  If  on  the  other  hand,  the  spectra  are  scaled  on  the 
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inertial  subrange  values  (c.f.  Kaimal  et  al.  1972)  it  is  clear  that  the  low  frequency  be¬ 
haviour  exhibits  significant  (order  of  magnitude)  departure  from  the  universal  curves. 

In  Figure  14,  the  uw  cospectra  for  the  same  runs  are  shown.  Recalling  that  z  and  U 
are  similar  for  all  runs,  it  is  clear  that  there  is  generally  substantial  departure  from  the 
imiversal  cospectrum  established  earlier.  In  particular  at  low  frequencies,  those  below 
that  of  the  swell,  the  stress  is  greatly  reduced,  or  even  upward.  This  reduction  in  stress 
is  consistent  with  the  general  reduction  in  turbulent  energy  at  these  frequencies,  as 
noted  above.  Reider  and  Smith  (1998),  in  their  study  of  the  effects  of  swell  on  the  stress 
vector,  concluded  that  the  major  effects  of  swell  on  stress  were  contained  in  the  mid¬ 
range  of  frequencies,  roughly  0.06-0.16  Hz.  Our  results  do  not  support  this,  but  indicate 
an  effect  at  all  frequencies  lower  than  that  of  the  swell. 

The  behaviour  described  above  is  least  evident  with  run  87056  (designated  by  o  in 
Figs.  12-14).  This  run  differs  from  the  remaining  fast  swell  runs  in  several  ways:  The 
inverse  wave  age  is  higher  than  that  of  the  other  cases  (0.61  versus  0.44-0.48)  and  the 
swell  energy  is  lower  by  roughly  a  factor  of  two.  Also,  the  swell  peak  frequency  is  higher. 
Designating  the  peak  wavenumber  of  the  swell  by  kp,  exp(—2kpz)  =  0.009  at  12m  for 
run  87056,  versus  0.06  for  the  remaining  runs.  Given  that  the  influence  of  the  waves  on 
the  atmospheric  turbulence  spectra  is  expected  to  decrease  with  height  as  exp(—2kpz) 
(Makin  and  Mastenbroek,  1996),  each  of  these  factors  would  reduce  the  observed  effect. 

Finally,  we  present  data  typical  of  slower  swell  cases,  those  with  U/cp  «  0.7.  Due  to 
the  fetch  geometry  of  the  lake  (Fig.  2),  there  is  generally  little  variation  in  the  swell 
frequency.  Hence  slow  swells  are  those  coexisting  with  higher  velocity  winds.  Here 
the  swell  is  the  result  of  a  decrease  in  wind  speed  (from  13  ms”^  to  7ms“^)  and  a 
40®  shift  in  wind  direction  which  took  place  just  prior  to  the  first  run.  The  resulting 
directional  wave  spectra  are  bimodal,  with  dominant  swell  peaks.  The  velocity  spectra 
and  cospectra  axe  plotted  along  with  the  pure  wind  sea  curves  in  Fig.  15.  The  data 
are  described  in  Table  2c.  Although  there  is  some  scatter,  in  general  the  data  are  in 
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agreement  with  the  universal  curves,  with  the  exception  of  Suw  for  run  87117.  Here, 
although  the  shape  of  the  curve  appears  consistent  with  the  universal  curve,  there  are 
significant  positive  stress  contributions  at  low  frequencies.  These  contributions,  perhaps 
associated  with  a  20®  wind  shift,  which  occurred  during  the  run,  lead  to  a  decreased  u*, 
and  hence  the  positive  offset  of  the  curve.  In  general  though,  the  slow  swell  cases  are 
consistent  with  the  pure  wind  sea  data.  This  might  be  anticipated  from  Fig.  10  in  that 
there  is  no  distinction  in  drag  between  the  pure  wind  sea  and  following  swell  cases  at 
these  wind  speeds. 

6.  Discussion 

We  have  presented  evidence  that  swell  effects  can  be  felt  at  a  considerable  height 
above  the  surface.  Velocity  spectra  and  cospectra  measured  over  pure  wind  sea,  with 
waves  of  a  wide  range  of  development  {U/cp  w  0.8  -  4.5),  agreed  with  the  universal 
cinrves  of  Miyake  et  al.  (1970).  Furthermore  measured  wind  velocities  and  spectra  at 
heights  as  low  as  2  m  {Hs  =  1.4m)  showed  no  evidence  of  significant  wave-coherent 
components.  This  is  consistent  with  the  estimates  of  Janssen  (1989)  and  Makin  and 
Mastenbroek  (1996),  and  implies  that  turbulent  stress  dominates  the  total  stress  at 
typical  measurement  heights  over  pure  wind  seas.  Likewise,  turbulent  spectra  taken  over 
slower  swells  {U/cp  w  0.7)  agreed  with  the  universal  curves.  Hence  Monin-Obukhov 
similarity  theory  is  expected  to  perform  well  in  these  conditions.  In  contrast,  with  fast, 
strong  swells  present  {U/cp  «  0.5,  Ha-swell  >  Ha-windsea)  wave  effects  were  evident 
in  the  wind  velocity  spectra  at  heights  of  12  m,  with  Hs  «  1.2m.  Turbulence  levels 
were  significantly  reduced  over  pure  wind  sea  values,  and  significant  deviations  from  the 
tiniversal  spectral  shapes  were  found.  Consequently,  we  do  not  expect  MO  theory  to  be 
valid  in  these  conditions. 

Smedman  et  al.  (1994)  recorded  fluxes  and  profiles  from  aircraft  and  a  nearshore 
tower,  and  observed  several  events  in  which  a  light  onshore  wind  was  blowing  shortly 
after  the  passage  of  a  storm.  Although  they  had  no  supporting  wave  data,  it  was 
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assumed  that  conditions  were  similar  to  those  observed  here,  with  swell  propagating 
onshore.  Their  observed  profiles  indicated  nearly  constant  wind  speed  with  height, 
consistent  with  their  estimations  of  very  low  surface  stress  (Cj)  ^  0-5  X  10-®).  They 
interpreted  the  observations  as  ‘shearfree’  turbulence,  with  no  active  shear  production  at 
the  air-sea  boundary.  Instead,  pressure  transport  was  found  to  be  significant  in  carrying 
turbulent  energy  produced  at  higher  elevations  by  shear  down  to  the  surface  layer.  Such 
a  mechanism  would  explain  the  apparent  paradox  observed  in  the  present  data:  that 
the  high  turbulent  energy  levels  in  the  inertial  subrange  are  apparently  not  supported 
by  the  classical  cascade  of  turbulent  energy  from  longer  scales.  Although  such  high 
level  pressure  measurements  were  not  made  during  the  WAVES  campaign,  near  surface 
measurements  show  an  upward  pressure  transport  from  the  swell  to  the  atmosphere 
during  fast  following  swell  conditions  (Kahma  et  al.  1999).  Hence  pressure  transport  is 
likely  important  in  the  presence  of  swell. 

As  was  pointed  out  by  Smedman  et  al.  (1994),  is  no  longer  a  valid  scaling  variable 
when  the  stress  is  near  zero  (or  negative).  Given  that  measured  stress  is  near-zero  in 
conditions  of  strong  following  swell,  it  is  perhaps  not  surprising  that  velocity  spectra 
and  cospectra  do  not  follow  MO  scaling  under  these  conditions.  There  are  significant 
consequences  arising  from  this.  First  of  all,  the  use  of  MO  based  stability  relations  for 
transforming  measured  stress  values  to  neutral  10  m  values  should  be  re-examined  for 
use  over  the  ocean.  Secondly,  the  commonly  used  inertial  dissipation  (ID)  method  for 
estimating  fluxes  must  also  be  called  into  question  in  use  over  the  open  ocean.  The 
method  is  derived  from  the  turbulent  kinetic  energy  (TKE)  conservation  equation 

odU  (uj0 -H  O.eiToW)  dm  dwpfp 

- z - ‘--dT-sr  =  ® 

where  e  and  p  are  fluctuations  in  turbulent  kinetic  energy  and  pressure  respectively. 

The  terms  are,  respectively,  production  of  TKE  from  shear,  production  (or  loss)  due  to 


-18- 


buoyancy,  loss  due  to  dissipation,  transport  of  TKE  and  pressure  transport.  Dividing  (7) 
by  uIJkz,  we  arrive  at 

KZ 

<f>u-C - =  0  (8) 

where  <j>t  represents  the  combined  transport  terms.  Following  MO  theory,  and  <f>t 
are  expected  to  be  universal  functions  of  C,  although  it  is  often  assumed  that  <f>f  =  0 
(Large  and  Pond,  1981).  Then  from  a  known  dissipation  rate  (usually  estimated  from 
the  inertial  subrange  of  the  frequency  spectrum,  assuming  Taylor’s  hypothesis),  and 
known  dimensionless  functions  (e.g.  those  of  Hogstrom,  1996),  the  stress  can  be  readily 
determined. 

The  method  has  gained  considerable  popularity  because  it  relies  only  on  the  high 
frequency  velocity  components,  which  are  not  as  susceptible  to  platform  motion  or  flow 
distortion.  Consequently  the  method  is  often  used  from  ships  at  sea  (Large  and  Pond 
1981;  Yelland  and  Taylor  1996).  Studies  comparing  momentum  fluxes  determined  by 
the  eddy-correlation  (EC)  and  ID  methods  have  typically  reported  good  agreement, 
with  rms  differences  of  approximately  20%  (Large  and  Pond  1981;  Edson  et  al.  1991). 
Donelan  et  al.  (1997)  found  similar  agreement  for  their  pme  wind  sea  data,  but  found  a 
doubling  of  the  rms  difference  in  their  comparison  when  swells  were  present.  When  fast, 
strong  opposing  swells  were  present,  the  EC  fluxes  were  approximately  double  the  ID 
fluxes. 

However,  <f>u  can  not  be  assiuned  to  be  a  function  of  ^  alone  in  the  event  that  the 
assumptions  of  MO  theory  are  violated.  Hence  on  the  basis  of  this  work,  ID  theory  can 
not  be  expected  to  work  in  the  presence  of  fast  swells.  This  is  evident  in  the  fast  swell 
data  of  Table  2b:  Whereas  the  eddy  correlation  data  exhibit  very  low  stress,  the  inertial 
dissipation  stress  values  differ  little  from  the  pure  wind  sea  values  (e.g.  nm  87176).  As 
seen  in  Figure  12,  turbulence  levels  in  the  inertial  subrange  are  similar  for  pure  wind 
sea  and  swell  runs  of  similar  wind  speed;  turbulence  levels  at  the  lower  frequencies 
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are  however  significantly  different.  When  spectra  do  not  follow  universal  scaling,  high 
frequency  behaviour  is  no  longer  a  good  surrogate  for  the  overall  spectrum. 

It  is  interesting  to  note  that  of  the  previous  work  on  the  effects  of  swell  on  momentmn 
fiux,  the  studies  reporting  no  significant  effect  were  based  on  inertial  dissipation  data 
(Dobson  et  al.  1994;  Yelland  and  Taylor  1996).  On  the  other  hand  an  effect  was  reported 
in  the  eddy  correlation  studies.  In  the  context  of  the  above,  it  is  not  surprising  that  this 
shovild  be  the  case.  It  is  now  clear  that  the  inertial  dissipation  method  itself  is  invalid  in 
these  conditions. 

7.  Conclusions 

1)  The  momentum  fiiix  can  be  strongly  affected,  even  to  the  extent  of  changing 
direction,  by  the  presence  of  swell.  Hence,  any  relationship  between  Cj)  and  U  must  take 
this  variable  into  account. 

2)  In  pure  wind  sea  conditions  spectra  follow  universal  scaling,  hence  Monin- Obukhov 
similarity  theory  would  appear  to  hold  [at  typical  measurement  heights].  This  is 
consistent  with  WBL  height  estimates  of  Janssen  (1969)  and  Makin  and  Mastenbroek 
(1996),  and  measurements  of  Edson  and  Fairall  (1998). 

3)  In  the  presence  of  fast  swells,  turbulent  velocity  spectra  and  cospectra  do  not  follow 
Monin-Obukhov  scaling.  In  particular,  turbulence  is  dramatically  reduced  at  frequencies 
about  and  lower  than  that  of  the  following  swell. 

4)  The  inertial  dissipation  method  of  estimating  fluxes  is  not  expected  to  be  valid  in  the 
presence  of  fast  swells  travelling  either  with  (this  paper)  or  against  (Donelan  et  al  1997) 
the  wind. 

5)  The  standard  practice  of  transforming  meteorological  data  to  neutral  conditions  is 
invalid  in  the  presence  of  fast  swells. 

In  midlatitudes,  such  fast  swells  are  often  associated  with  frontal  passages,  or  may 
occiur  following  a  storm.  In  the  tropics,  where  the  winds  are  weaker,  such  conditions 
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often  occtir  in  advance  of,  or  following,  the  passage  of  tropical  storms.  Hence  the  results 
have  particular  significance  in  accurate  modelling  of  extreme  events. 

For  many  years,  campaigns  to  determine  the  drag  coefficient  over  the  sea  have 
been  carried  out.  Yet  the  experimental  scatter  in  the  data,  both  within  and  between 
experiments,  remains  large,  especially  at  low  wind  speeds.  The  results  here  suggest  that 
this  is  due,  at  least  in  part,  to  the  presence  of  swells.  Other  conditions,  such  as  rapidly 
turning  winds,  are  also  likely  factors  (Rieder  1997).  In  order  for  bulk  formulations  to 
successfully  reflect  these  factors,  Monin-Obukhov  similarity  theory  must  be  extended  to 
include  additional  variables.  This  study  suggests  that  the  swell  related  variables  U/cp, 
kpZ  and  fTa-swell/ffa-windsea  are  good  candidates,  although  clearly  further  work  is 
required  to  quantify  the  effects. 
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FIG.  8.  Representative  wind  sea  data  in  the  universal  scaling  of  Miyake  et  al.  (1970). 
The  panels  show  spectra  of  vertical  velocity  (a)  and  cospectra  between  the  vertical 
and  horizontal  in-line  velocity  (b).  The  data  represented  by  o  and  x  axe  from  the  12  m 
bivane  with  wind  speed  U  —  9.6  and  4.7ms“^respectively.  The  *  data  are  from  a  hotfilm 
at  2.1m,  during  the  same  run  as  the  o  data  -  see  Table  2a.  The  shaded  area  denotes  ±1 
standard  error  from  the  mean  curve  established  from  all  wind  sea  data. 

FIG.  9.  Stress  offwind  angle  with  respect  to  the  mean  wind  direction  vs.  12  m  wind 
speed.  Pure  wind  sea  (o);  mixed  seas  (x). 

FIG.  10.  12  m  non-neutral  data  from  WAVES  (pme  wind  sea,  •,  and  following  swell, 
x)  and  SWADE  (Donelan  et  al.  1997:  pure  wind  sea,  *,  and  counter  swell,  o). 

FIG.  11.  Spectra  from  run  87146,  with  light  winds  and  fast  following  swell.  Curves  are 
downwind  (xlO)  and  vertical  velocity  spectra  from  bivane  at  11.9m  and  svirface  elevation 
spectrum.  U  —  3.9ms“^,  Hs  =  1.2m,  and  U/cp  =  0.44. 

FIG.  12.  Spectra  of  vertical  (a)  and  horizontal  in-line  (b)  velocity  during  swell  runs. 
The  data,  described  in  Table  2b,  are  from  the  wind  speed  range  C/"  «  4  —  5ms~^.  For  the 
o  data,  j7/cp  =  0.61;  for  the  other  data,  U/cp  <  0.5.  The  solid  line  is  the  pure  wind  sea 
universal  curve,  rendered  dimensional  using  the  data  of  pure  wind  sea  run  87176. 

FIG.  13.  Spectra  of  vertical  (a)  and  horizontal  in-line  (b)  velocity  from  swell  runs, 
plotted  in  the  imiversal  scaling  of  Miyake  et  al.  (1970).  The  data,  described  in  Table  2b, 
are  from  the  wind  speed  range  U  w  4  — 5ms~^.  For  the  o  data,  U/cp  =  0.61;  for  the  other 
data,  U fcp  <  0.5.  The  shaded  areas  represent  ±  one  standard  error  about  universal 
curves  established  from  the  pure  wind  sea  data. 

FIG.  14.  Bivane  velocity  uw  cospectra  from  swell  runs.  The  data,  described  in  Table 
2b,  axe  from  the  wind  speed  range  U  Pi  A  —  5ms~^.  For  the  o  data,  U/cp  =  0.61;  for 
the  other  data,  U/cp  <  0.5.  The  solid  line  is  the  pure  wind  sea  universal  curve,  rendered 
dimensional  using  the  data  of  pure  wind  sea  run  87176. 
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FIG.  15.  Bivane  velocity  spectra  (vertical  (a)  and  horizontal  inline  (b))  and  cospectra 
(c)  from  slow  swell  runs,  in  the  universal  scaling  of  Miyake  et  al  (1970).  The  data  are 
described  in  Table  2c.  The  shaded  areas  represent  ±  one  standard  error  about  the  pure 
wind  sea  universal  curves. 
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